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Abstract: In the planar Af = 4 supersymmetric Yang-Mills theory, the conformal sym- 
metry constrains multi-loop n-edged Wilson loops to be basically given in terms of the 
one-loop n-edged Wilson loop, augmented, for n > 6, by a function of conformally invari- 
ant cross ratios. We identify a class of kinematics for which the Wilson loop exhibits exact 
Regge factorisation and which leave invariant the analytic form of the multi-loop n-edged 
Wilson loop. In those kinematics, the analytic result for the Wilson loop is the same as 
in general kinematics, although the computation is remarkably simplified with respect to 
general kinematics. Using the simplest of those kinematics, we have performed the first 
analytic computation of the two-loop six-edged Wilson loop in general kinematics. 
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1. Introduction 



In the planar N = 4 supersymmetric Yang-Mills (SYM) theory, Anastasiou, Bern, Dixon 
and Kosower (ABDK) [1] proposed an iterative structure for the colour-stripped two-loop 
scattering amplitude with an arbitrary number n of external legs in a maximally-helicity 
violating (MHV) configuration. Writing at any loop order L, the amplitude as the 

tree-level amplitude, Mn \ which depends on the helicity configuration, times a scalar 
function, m[ i ', 



Thus the two-loop amplitude is determined in terms of the one-loop MHV amplitude m n (e) 
evaluated through to 0(e 2 ) in the dimensional-regularisation parameter e = (4 — d)/2, the 
constant = -Cf/ 2 > an d the function f( 2 \e) = -(2 - C3^ - C^ 2 , with & = CW and ((z) 
the Riemann zeta function. 

Subsequently, Bern, Dixon and one of the present authors (BDS) proposed an all-loop 
resummation formula [2] for the colour-stripped n-point MHV amplitude, which implies a 
tower of iteration formulae, allowing one to determine the ra-point amplitude at a given 
number of loops in terms of amplitudes with fewer loops, evaluated to higher orders of 
e. BDS checked that the ansatz is correct for the three-loop four-point amplitude, by 
evaluating analytically mf\e) through to finite terms, as well as m^\e) through to 0(e 2 ) 
and m'p {e) through to C(e 4 ). The ansatz has been proven to be correct also for the two- 

(2) 

loop five-point amplitude [3, 4], for which (e) has been computed numerically through 
to finite terms, as well as m^\e) through to 0{e 2 ). 

Using the AdS/CFT correspondence, Alday and Maldacena showed that in the strong- 
coupling limit the ansatz must break down for amplitudes with a large number of legs [5]. 
At weak coupling, the computation of the two-loop six-edged Wilson loop [6] led to conclude 
that either the ansatz or the duality relation between amplitudes and Wilson loops were 
to break down for two-loop six-point amplitudes. Likewise, there were hints of a failure 
of the ansatz from the six-point amplitude analysed in the multi-Regge kinematics in a 
Minkowski region [7, 8, 9]. The accumulating evidence against the ansatz provoked the 
numerical calculation of rn^\e) through to finite terms and of m^\e) through to 0(e 2 ), 

by which the ansatz was demonstrated to fail [10], and where it was shown that the finite 

(2) 

pieces of the parity-even part of m 6 (e) are incorrectly determined by the ansatz (although 

(2) 

the parity-odd part of m & (e) does satisfy the ansatz [11]). In particular, it was shown 
numerically that the two-loop remainder function, defined as the difference between the 
two- loop amplitude and the ansatz for it, 




(1.1) 



™P (e) = - [m« (e)J + (e) m« (2e) + + 0(e) . 



(1.2) 




1 r 



(1.3) 
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is different from zero for n = 6, where Rn is a function of the kinematical parameters of 
the n-point amplitude, but a constant with respect to e. However, the analytic form of 

(2) 

R 6 was not computed. 

In the strong-coupling limit, Alday and Maldacena [12] showed that planar scattering 
amplitudes exponentiate like in the ansatz, and suggested that the vacuum expectation 
value of the n-edged Wilson loop could be related to the n-point amplitude in TV = 4 
SYM. At weak coupling, the agreement between the light-like Wilson loop and the (parity- 
even part of the) MHV amplitude has been verified for the one-loop four-edged [13] and 
n-edged [14] Wilson loops, and for the two-loop four-edged [15], five-edged [16] and six- 
edged [6, 17] Wilson loops. 

Furthermore, it was shown that the L-loop light-like Wilson loop exhibits a conformal 
symmetry, and that the solution of the Ward identity for a special conformal boost is 
the ansatz, augmented, for n > 6, by a function R^ly L of conformally invariant cross- 
ratios [16]. Because of the duality between Wilson loops and amplitudes at one and two 
loops, the function of the conformally invariant cross-ratios, R^lvi' can ^ e identified as 
the remainder function of Eq. (1.3). 

In Refs. [17, 18], the two-loop n-edged Wilson loop has been given in terms of Feynman- 
parameter-like integrals. Furthermore, in Ref. [18] a numerical algorithm has been set up, 
which is valid for the two-loop n-edged Wilson loop and by which the seven-edged and eight- 
edged Wilson loops have been computed (although the corresponding MHV amplitudes are 
not known 1 ). Thus, also the remainder functions R i f\ VL and R^ WL of the Wilson loops 
are known numerically, and the numerical evidence [18] confirms that they are functions of 
conformally invariant cross-ratios only. However, their analytic form is in general unknown. 

In fact, before this letter the remainder function had been computed analytically only 
in the strong coupling regime: for the two-loop six-edged Wilson loop in the particular 
kinematic configuration in which the cross ratios coincide [20] , and otherwise in a particular 
kinematic set-up for which only 2n-edged polygons are allowed [21]. In that set-up, the 
simplest non-trivial remainder function is the one of the two- loop eight-edged Wilson loop. 
In the weak coupling regime, and in the same kinematics, R^ WL has been computed 
numerically in Ref. [22], where a numerical evidence has been found for a linear relation 
between the remainder functions at weak and strong couplings. 

In this letter, we give a brief account of the first analytic computation at weak coupling 
of the two-loop six-edged Wilson loop in general kinematics. The computation has been 
done in the Euclidean region in D = 4 — 2e dimensions, where the result is real. 

In Sec. 2, we write the two-loop Wilson loop in terms of the one-loop Wilson loop 

(2) (2) 

plus a remainder function R y n ' WL . For n = 6, R 6 WL is a function of the three conformally 
invariant cross ratios, ni,U2,U3. It has been computed numerically in Refs. [10, 17, 18] in 
arbitrary kinematics, constrained only by momentum conservation. However, it suffices to 
compute R^ WL in any kinematical limit which does not modify the analytic dependence of 

(2) 

Rq y VL on ui, U2, us. To that end, we recall that the L-loop four-edged Wilson loop is not 

lr The parity-even part of the two- loop n-point MHV amplitude has been given in terms of integral 
functions, yet to be evaluated [19]. 
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modified by the Regge limit [13]. In such a limit w 4 undergoes an exact Regge factori- 
sation. So is the case of the five-edged Wilson loop in the multi- Regge kinematics [7, 23]. 
However, the six-edged Wilson loop is modified by the multi-Regge kinematics: the three 
conformally invariant cross ratios are not invariant in such a limit [23]. Less constraining 
Regge limits have been analysed in Ref. [24]. The simplest of those limits to feature an 
exact Regge factorisation of is the quasi-multi-Regge kinematics (QMRK) of a pair 
along the ladder [25, 26]. 

In Sec. 3, we recall the QMRK of a pair along the ladder for the six-edged Wilson loop, 
and we show that the QMRK of three-of-a-kind along the ladder [27] for the seven-edged 
Wilson loop, the QMRK of four-of-a-kind along the ladder [28] for the eight-edged Wilson 
loop, and in general the QMRK of a cluster of (n — 4)-of-a-kind along the ladder for the 
n-edged Wilson loop do not modify the analytic dependence of vffl on the conformally 
invariant cross ratios. That is, this class of kinematics exhibits an exact Regge factorisation 
of w^h \ Thus, the result for in these kinematics is the same as the result in general 
kinematics, although the computation is remarkably simplified with respect to the same 
computation in general kinematics. Finally, we note that although in Sec. 4 we apply the 
analysis of Sec. 3 to the computation of the six-edged two-loop Wilson loop, nothing of 
what we consider in Sec. 3 is specific to two loops: The analysis of Sec. 3 is valid for any 
number of loops. 

In Sec. 4, we brief on how the Feynman-parameter-like integrals of the two-loop six- 
edged Wilson loop have been computed in the QMRK of a pair along the ladder, and on the 
type of functions which appear in the final result. Because of the exact Regge factorisation, 
the ensuing remainder function is valid in general kinematics. It can be expressed as a linear 
combination of multiple polylogarithms of uniform transcendental weight four. However, 
the result is far too long to be reported in this letter. We present it in an electronic form at 
www . arxiv . org where a text file containing the Mathematica expression for the remainder 
function is provided. 

2. The two- loop Wilson loop 

The Wilson loop is defined through the path-ordered exponential, 



computed on a closed contour C n . In what follows, the closed contour is a light-like n-edged 
polygonal contour [12]. The contour is such that labelling the n vertices of the polygon as 
xi, . . . ,x n , the distance between any two contiguous vertices, i.e., the length of the edge 
in between, is given by the momentum of a particle in the corresponding colour-ordered 
scattering amplitude, 



with i = l,...,n. Because the n momenta add up to zero, Y^i=iPi = 0, the n-edged 
contour closes, provided we make the identification x\ = x n+ \. 




(2.1) 



Pi — %i Xi-\-\ , 



(2.2) 
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In the weak-coupling limit, the Wilson loop can be computed as an expansion in 
the coupling. The expansion of Eq. (2.1) is done through the non-abelian exponentiation 
theorem [29, 30], which gives the vacuum expectation value of the Wilson loop as an 
exponential, 

oo oo 

[W[C n ]) = 1 + a L W^ = exp «M L) , (2-3) 

L=l L=l 

where the coupling is defined as 

For the first two loop orders, one obtains 

v>P=WP, wV = W<?)-±(wU) 2 • (2.5) 

The one-loop coefficient was evaluated in Refs. [13, 14], where it was given in terms 
of the one-loop n-point MHV amplitude, 

«£> = = mi n ] - 4 + w . ( 2 - 6 ) 

where the amplitude is a sum of one- loop two-mass- easy box functions [31], 

m (i) = J2F 2 ™(p ,q,P,Q), (2.7) 

where p and q are two external momenta corresponding to two opposite massless legs, 
while the two remaining legs P and Q are massive. The two-loop coefficient has been 
computed analytically for n = 4 [15] and n = 5 [16] and numerically for n = 6 [17] and 
n = 7,8 [18]. 

In Ref. [16] it was established that the Wilson loop fulfils a special conformal Ward 
identity, whose solution is the ansatz plus, for n > 6, an arbitrary function of the confor- 
mally invariant cross-ratios, defined in Eq. (2.11). Thus, the two-loop coefficient Wn^ can 
be written as 

^i 2) (e) = (26) + C$ L + R% L + 0(e) , (2.8) 

where the constant is the same as in Eq. (1.2), C$ L = = -Cl/2, and the function 
f$ L (e) is [15, 18, 32] 

Zv?L(e) = -C2 + 7C 3 e-5C4e 2 . (2.9) 

With the two-loop coefficient Wn^ given by Eqs. (2.8) and (2.9) and the two-loop MHV 
amplitude given by Eqs. (1.2) and (1.3), the duality between Wilson loops and amplitudes 
is expressed by the equality of their remainder functions [17, 18], 

R ( nWL = R { n ] ■ (2-10) 



n,WL 

Defining the conformally invariant cross ratios as, 



x ij x i+ij+i 
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for n = 6, they are [16], 

„, _ „, _ x 13 x 46 „ _ „, _ x 15 x 24 „, _ „, _ x 26 x 35 /o -i o\ 

^36 — «1 — —o o , «14 — ^2 — —9 o- , U25 — U3 — —9 9- , I 2 - 12 / 

x 36 x 41 x 14 x 25 J/ 25 x 36 

where xfj = (x{ — Xj) 2 , and using Eq. (2.2) one sees that xf i+2 = Sj,i+i and x 2 i+3 = 
Si,i+i,i+2, where the labels are understood to be modulo 6. 

3. The quasi-multi-Regge kinematics of a cluster along the ladder 

(2) 

As we remarked in the Introduction, it suffices to compute R 6 WL in any kinematical limit 

(2) 

which does not modify the analytic dependence of R\ ! WL on ui, 1*2, 113. The simplest of 

those limits to feature an exact Regge factorisation of and in fact in general of 
is the QMRK of a pair along the ladder [25, 26]. In those kinematics, the outgoing gluons 
are strongly ordered in rapidity, except for a central pair of gluons along the ladder, while 
their transverse momenta are all of the same size. In the physical region, defining 1 and 2 
as the incoming gluons and 3, 4, 5, 6 as the outgoing gluons, the ordering can be chosen as 

2/3 > 2/4 ^ 2/5 > ^6; |P3±| - \P4±\ - \P5±\ - \P6±\ , (3-1) 

where the particle momentum p is parametrised in terms of the rapidity y and the azimuthal 
angle (/>, p = (\p± \ coshy, \p± \ cos</>, \p± \ sm<ft, \p± \ sinhy). We shall work in the Euclidean 
region, where the Wilson loop is real. There the Mandelstam invariants are taken as all 
negative, and in the QMRK of a pair along the ladder they are ordered as follows, 

Sl2 > -«34, -S56, -S345, -S123 > -S23, -S45, S61, -S234 • (3.2) 

Introducing a parameter A <C 1, the hierarchy above is equivalent to the rescaling 

{^34,556, 5123,5345 } = O(X) , {s 2 3 , ^45, «61 , «234} = 0(\ 2 ) . (3.3) 

It is easy to see that in this limit the three conformally invariant cross-ratios (2.12) do not 
take trivial limiting values [24], 

Ui _^ U QMRK _ £45 = 

1 (Pi +Pt)(Pl +Pb) 

^ ^ u QMRK = \P3±\ 2 pjp e = m j 

(\P3± +P4_l| 2 +ptPi)ipt +Pb )P% 
|2„+„ 



QMRK = \P6±FptP: = 0{l) . 

pt(P4 +P 5 )(IP3_L +P4±| 2 +P^>4 ) 



A similar analysis can be carried through for the seven-edged Wilson loop, w\ L \ We 
have verified that the simplest limit to feature an exact Regge factorisation is the QMRK 
of three-of-a-kind along the ladder [27]. In the physical region, the outgoing gluons are 
strongly ordered in rapidity, except for a cluster of three along the ladder, 

2/3 > 2/4 - 2/5 - 2/6 > VT, \P3±\ - \P4±\ - \P5±\ - \Pe±\ - \P7±\ ■ (3.5) 
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(3-7) 



In the Euclidean region, the Mandelstam invariants are ordered as follows, 

—S12 > -«123, -^345, -^567, -^712, -S34, -«67 > .„ „. 

(3.6) 

> -S23, -S45, —S56, —$71, — «234, — S456; — «671 • 

Through a parameter A <C 1, the hierarchy above is equivalent to the rescaling 

{S123, S345, S 567 , S712, S34, S 67 } = 0(A) 
{>23, S45, «56, S71, S 2 34, S 4 56, S 6n } = 0{\ 2 ) . 

Using Eq. (2.11), and the fact that a;|- +1 = Sj...j, it is easy to see that the seven cross ratios 
of the seven-edged Wilson loop do not take trivial limiting values under the rescaling (3.7), 

{uu,U25,U 3 q,U4 7 ,U 5 i,Uq 2 ,U 73 } = 0(1) . (3.8) 

Thus, the dependence of on the seven cross ratios is not modified by the QMRK of 
three-of-a-kind along the ladder (3.6), and hence undergoes an exact Regge factorisa- 
tion in this limit. 

The same pattern unfolds for the eight-edged Wilson loop, . The simplest limit to 
feature an exact Regge factorisation is the QMRK of four-of-a-kind along the ladder [28]. 
In the physical region, the outgoing gluons are strongly ordered in rapidity, except for a 
cluster of four along the ladder, 

2/3 > 2/4 ^ 2/5 ^ 2/6 ^ 2/7 > V8\ \P3± \ - \P4±\ - \P5±\ — \P6±\ - \P7±\ - \P8± \ ■ (3.9) 

In the Euclidean region, the Mandelstam invariants are ordered as follows, 

— S\2 3> -S1234, — •S3456, —$123, -S345, -S678; ~ «812, ~ ^34, ~ «78 > 

~> -S2345; —•54567, ~S234, -S456, -«567 5 ~ $781, ~S23, ~S45, ~$56, -$67, ~ «81- 

Through the parameter A <C 1, the hierarchy above corresponds to the rescaling 

{^1234, S3456, «123, ^345, S678, s 812, S34, Sfg} = 0(A) , 
{•S2345, S4567, ^234, «456; ^567, ■5781, s 23, S4b, S56, S67, S8l} = 0{\ 2 ) . 



(3.10) 



(3.11) 



It is easy to check that all the twelve cross ratios of the eight-edged Wilson loop do not 
take trivial limiting values under the rescaling (3.11), 

{UU,U25,U36,U4 7 ,U 5 S,U61,U 7 2,U 8 3,U1 5 ,U26,U3 7 ,U4 8 } = 0(1) . (3-12) 

Thus, the dependence of on the twelve cross ratios is not modified by the QMRK of 
four-of-a-kind along the ladder (3.10), and hence undergoes an exact Regge factorisa- 
tion in this limit. 

The pattern above generalises to the n-edged Wilson loop, We illustrate briefly 

how the QMRK of a cluster of (n — 4)-of-a-kind along the ladder features the exact Regge 
factorisation of . In the physical region, the outgoing gluons are strongly ordered in 
rapidity, except for a cluster of (n — 4) along the ladder, 

2/3 » 2/4 ^ ■ ■ ■ ^ Vn-\ » Vn, \P3±\ - ■ ■ ■ - \Pn±\ ■ (3-13) 



-6- 



In order to display the exact Regge factorisation of we deal with the cases of an 

even, n = 2r, and an odd, n = 2r + 1, number of edges separately. Through a parameter 
A <C 1, the hierarchy of the Mandelstam invariants implied by Eq. (3.13) can be rendered 
by requiring that S12 = 0{1) and by the rescaling 

s 12 ■ = 0(A), 3 < j < r, 

1 v ; , (3.14) 
S34...J =C(A), 4<j <r + 2, 

for any n and r > 3. In addition, 

s i2r =0(X), r + 2<7<2r-l, r>3, 

J V 7 ~ ~ (3.15) 

Sj...i2 = 0(A) , r + 4<j<2r, r>4, 

for n = 2r, with the labels understood to be modulo 2r, and 

Si...2r+i = O(A) , r + 2<j<2r, 

(3.16) 

83...12 = 0{\) , r + 4 < j < 2r + 1 , 

for n = 2r + 1 and r > 3, with the labels understood to be modulo 2r + 1. All other 
invariants rescale to be C(A 2 ). It is easy to check that for n = 6,7,8, Eqs. (3.14)-(3.16) 
reproduce the rescaling of Eqs. (3.3), (3.7) and (3.11). 

In order to compare with the rescaling of Eqs. (3.3), (3.7) and (3.11), it is convenient 
to classify the cross ratios as follows, 

Ulj = = , J = 4,...,r + 1, 

X lj X 2j+l Si...j-iS2~-j 

u 2j = X 4^ = S -^^, j = 5,...,r + 2 > 

x 2j x ij+i s 2 ...j-is 3 ...j ;:5.17) 



_ x 2r+lj+l x lj _ S 2r +l-jSl...j-l 

u 2r +ij — — 2 — ' J — o, . . . , r , 

X 2r+lj X lj+l S2r+l-j-lSl—j 

for n = 2r + 1, and 

», 4^ .y 1 r-l. 

X lj X 2j+l S\...j-\S 2 ...j 



2 2 

c rj+l x r+lj _ S r ...jS r +i...j-i 



X rj X r+lj+l S r ...j-iS r+ i. 



j = r + 3, . . . , 2r , 



x r+lj+l x r+2j S r+ i...jS r+2 ...j-i . (3.18) 

u r+ij = 2 2 = — ' j = r + 4, . . . , 2r , 

X r+lj X r+2j+l s r+l---j-lS>r+2-j 



_ X lrj+l X lj _ S 2r ... jSl ...^i 
a 2rj — 2 2 — ' — o, . . . , / 

X 2rj X lj+l S 2r ...j-\S\...j 
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for n = 2r, where we use the fact that u^+i = = 0. Without further imposing the 

Gram-determinant constraints, the counting yields n(n — 5)/2 conformally invariant cross 
ratios, in agreement with Ref. [18]. It is easy to check that for n = 6,7,8, Eqs. (3.17) 
and (3.18) generate the cross ratios of Eqs. (2.12), (3.8) and (3.12). Then by inspection 
one can see that the cross ratios (3.17) and (3.18) do not take limiting values under the 
rescaling of Eqs. (3.3), (3.7) and (3.11). Thus, the dependence of on the cross ratios 
(3.17) and (3.18) is not modified by the QMRK of a cluster of (n — 4)-of-a-kind along the 
ladder (3.13), and hence undergoes an exact Regge factorisation in this limit. 

Finally, we note that the exact Regge factorisation of the n-edged Wilson loops, with 
n = 4, 5, may be dealt with as a degenerate case of the QMRK of a cluster of (n — 4)-of- 
a-kind along the ladder. Namely, for n = 4 one obtains the QMRK of a cluster of zero 
particles along the ladder, which is the standard Regge limit [13], and for n = 5 the QMRK 
of a cluster of one particle along the ladder, which is the multi- Regge kinematics of five 
particles [7, 23]. 

4. The two-loop six-edged Wilson loop 

In Ref. [18] an expression for a generic two-loop n-edged Wilson loop as a sum of Euler- 
type integrals was presented, similar to Feynman parameter integrals appearing in the 
computation of Feynman integrals. In Sec. 3 we argued that the Wilson loops are Regge 
exact in the QMRK where (n — 4) gluons are emitted along the ladder. Hence, it follows 
that it is sufficient to compute the individual integrals in the QMRK to obtain the Wilson 
loop in general kinematics. 

In the present work we concentrate exclusively on the two-loop six-edged Wilson loop, 

(2) 

, which is the first case where the remainder function is non zero. Hence, an analytic 

computation of is equivalent to an analytic computation of the two-loop six-point 

(2) 

remainder function RyJ L 6 - We start from the parametric representations for the diagram- 
matic contributions to the two-loop n-edged Wilson loop derived in Ref. [18], and we derive 
appropriate Mellin-Barnes representations for them using the standard formula, 



where the contour is chosen such as to separate the poles in T(. . . — z) from the poles in 
T(. . . + z). Note that in our case A is in general an integer plus an off-set corresponding 
to the dimensional regulator e. In order to resolve the singularity structures in e, we apply 
the strategy based on the Mellin-Barnes representation and given in Refs. [33, 34, 35, 36]. 
To this effect, we apply the codes MB [37] and MBresolve [38] and obtain a set of Mellin- 
Barnes integrals that can be safely expanded in e under the integration sign. Then we 
proceed to take the QMRK limit defined by Eq. (3.2) by applying MBasymptotics [39] to 
extract the leading QMRK behavior of each Mellin-Barnes integral, and barnesroutines 
[40] to perform integrations that can be done by corollaries of Barnes lemmas. To be more 
explicit, in Ref. [18] the six-edged Wilson loop was expressed as, 




(4.1) 




(4.2) 



z 



8 



where fi(pk) denote the parametric integrals of Ref. [18] depending on the external mo- 
menta pk ■ The prefactor C is defined by 



C = 2aV e [r(l + e)eT B£ ] 2 . (4.3) 



and the scale /u 2 is given in terms of the Wilson loop scale, = Tre yEe fj 2 . The Regge 
exactness of enables us to write 

w^=Cj2t\Pk), (4-4) 

i 

where f^\pk) denotes the leading asymptotic behavior of fi(pk) in the QMRK limit defined 
by the scaling (3.3), 

h(Pk) = fi 1 \pk) + 0(\). (4.5) 

In Section 3 we considered the QMRK limit where gluons 1 and 2 are incoming, i.e., where 
s±2 is the largest invariant. Of course there are five additional ways in which we could have 
defined the limit, corresponding to the cyclic permutations of the external gluons, e.g., we 
could have considered the QMRK limit defined by the scaling, 

{S45,S61,S234,S123} = 0(A) , {s 34 , S 56 , S 12 , S 345 } = 0{\ 2 ) , (4.6) 

where A < 1 and S23 = 0(1). Note that this limit is incompatible with the limit (3.3), i.e., 
terms that are 0(A) in one limit can be large in another limit, and vice-versa. However, 
the Regge exactness of the Wilson loop allows us to iterate this procedure and to repeat 
the previous argument starting from Eq. (4.4) and to take the limit (4.6). Then we arrive 
at 

i 

where f\ 2 \pk) is the leading asymptotic behavior of f\ l \pk) in the limit (4.6), 

ri 1 \p k ) = fP(Pk) + 0(\). (4.8) 

It is straightforward to see how this procedure iterates for the remaining four cyclic permu- 
tations of Eq. (3.3). Finally, we arrive at a set of multiple Mellin-Barnes integrals f^ 6 \pk) 
of a much simpler type than the original ones. After applying our procedure, all integrals 
are at most threefold and all of them are explicitly dependent on the conformal cross ratios 
only, because the cross ratios are the only combination of invariants that are invariant 
under all the cyclic permutations of Eq. (3.3). However, note that the coefficients of the 
integrals do not depend only on the conformal cross ratios, but also on logarithms of Man- 
delstam invariants, which arise when expanding the MB integrals in a QMRK limit of the 
type defined in Eq. (3.3), as can be easily seen from the following example, 

-I p- l/2+ioo / \ z 

— \ dz T(-z) T(z) AM = In -H + 0(A) . (4.9) 
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Finally, we checked numerically that the sum of the Mellin-Barnes integrals in the QMRK, 
fi \Pk), is equal to the sum of all the original parametric integrals, the latter being eval- 
uated numerically using FIESTA [41, 42]. 

The resulting Mellin-Barnes integrals are then evaluated by directly closing contours 
and summing up residues or by exchanging a Mellin-Barnes integration with an integral of 
Euler type. The infinite sums which appear in the intermediate steps of the computation 
are typically generalised harmonic sums [43, 44] as well as multiple binomial sums [45, 46]. 
The convergence of the series requires the conformal cross ratios to be less than 1, and 
in the following we concentrate on this kinematic region, within the Euclidean region. 
Details on the explicit computation of the integrals will be presented in a forthcoming 
publication [47]. Here it suffices to say that, except for the contribution coming from 
the hard diagram with six light-like edges, all the integrals can be expressed in terms of 
harmonic polylogarithms [48] in one conformal cross ratio. In turn, the six-edged hard 
diagram constitutes the bulk of the final result, and can be written as a linear combination 
of Goncharov's multiple polylogarithms [49], whose arguments are functions of conformal 
cross ratios. These polylogarithms are defined by the iterated integration, 



where we define w = (a,w'), and for z = 1 they are manifestly real, if all the elements 
in the weight vector w are either greater than 1 or negative. The number of elements of 
w is called the (transcendental) weight of G(w;z). The polylogarithms we obtain can be 
divided into several classes, corresponding to the elements Wi of the weight vector, 

1. Wi = 1/Uj, 1/(1 - Uj), (1 - Uj)/(1 - Uj - u k ). 

It is easy to see that in this case > 1 or Wi < 0, for < Ui, Uj < 1. 

2. Wi = l/(ui + Uj). 

In this case Wi could be smaller than 1, i.e., the polylogarithms can develop an 
imaginary part. However, we checked numerically that the imaginary parts cancel in 
the final answer. 

3. Wi = l/vf£J, where we define 




(4.10) 



,(±) 
■jkl 



1-Uj -u k + ui±J (uj + u k - ui - l) 2 - 4 (1 - uj) (1 - u k ) ui 



u 



2(1- Uj) ui 




(4.11) 



v 



jkl 



2(1- Uj) u k 



A comment is in order about the square roots in Eq. (4.11). It turns out that the square 
roots become complex for certain values of the conformal cross ratios inside the unit cube, 
but they always come in pairs such that the sum of the two contributions is real. To 



emphasize this property, we introduce the following notation, 

uijk, ...; z) = G ug^, ...;*)+ G , u$,. . . ; zj , 
H(w; 1/uijk) = H (w; + H (w; l/u{ 



(4.12) 



and similarly for u,-^. 

After having computed the contributions from the individual integrals, we can easily 

(2) 

extract the remainder function R\y L 6 by subtracting the contribution from our computa- 

(2) 

tion, thus obtaining the first fully analytic representation of RyJ L 6 in the Euclidean region. 
Note that although we performed the computation in the QMRK of a pair along the ladder, 
as introduced in Sec. 3, the Regge exactness of the Wilson loop in this limit ensures that 
our expression is valid in general kinematics. The final result for the remainder function 
can be expressed as a linear combination of multiple polylogarithms of uniform transcen- 
dental weight four 2 . Because the result is rather lengthy, we present it in an electronic 
form at www.arxiv.org where a text file containing the Mathematica expression for the 
remainder function is provided. 

We have checked numerically that our result is completely symmetric in its arguments. 
Furthermore, we have checked analytically that the expression satisfies the constraints 
imposed by the multi-Regge and the collinear limits. Note that the vanishing in these limits 

(2) 

is non trivial, since the expression of RyJ L 6 in general kinematics involves polylogarithms 
whose arguments are ratios of cross ratios, which can be 0{1) in the limit. However, all 
those contributions exactly cancel when approaching the limit. Finally, we have checked 
numerically at several points that our results agree with the numerical results of Ref. [18]. 
In the particular case where all three conformal cross ratios are equal, we find that, 



R WL 6 ( 1 > !> !) = -?« ~ -a-''U!>Bl... , 

36 4 (4-13) 

J™ = ~Tu ~ -°- 67645 - ' 

in agreement, within numerical errors, with the values quoted in Ref. [18]. Similarly, the 
asymptotic behavior for u — > is given by, 



lim R$ L>6 (u, u,u) ^Uru ' ' ( ' ) 



Further results for the special case where all three conformal ratios are equal are summa- 
rized in Fig. 1. Note that even though our numerical evaluation is for the moment limited 
to < Ui < 1, we can still compute the asymptotic value when all conformal cross ratios 
are equal and large by expanding the Mellin-Barnes integrals around u = oo before taking 
residues. We find perfect agreement with the numerical value quoted in Ref. [18], which 
deviates slightly from the asymptotic value obtained from the analytic expression of the 
remainder function proposed in Ref. [20]. 



2 In the present version of the remainder function, we had to extend the definition of the transcendental 
weight to include the imaginary roots, which are present in Eq. (4.11). However, given that the remainder 
function is real, we cannot exclude that in a suitable basis those imaginary roots could disappear. 
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R(,(u, it. u) 




Figure 1: The remainder function i?g (u, u, u) for < u < 1. The points represent the numerical 
values given in Ref. [18]. 



We conclude this section by making some comments on the numerical evaluation of 
Goncharov's multiple polylogarithms. Up to weight two, Goncharov's polylogarithms can 
be expressed in terms of ordinary logarithms and dilogarithms, e.g. 3 , 

G(a;z) =ln(l-£) , 

0(o , M) . Li2 (^£)_ Li2 (_A_) +log(l _£ )lo /-^ ( " 5) 

However, our result involves polylogarithms up to weight four. We observe that in these 
cases the polylogarithms can be evaluated in a numerically stable and fast way, even for 
complex arguments, by writing them as an iterated integration of polylogarithms of weight 
two, 



G(a,b,c;z) = / - — — G(&,c;ii), 

G(a,b,c,d;z) = f dtl f -^—G(c,d;t 2 ), 
Jo h- a Jo t 2 -b 



(4.16) 



and using Mathematical native NIntegrate command to perform the integration, provided 
that the integrals converge. 



5. Conclusions 

In this letter we have identified a class of kinematics for which the multi-loop n-edged 
Wilson loop exhibits exact Regge factorisation and which leave invariant the analytic form 
of the Wilson loop. In those kinematics, the analytic result for the Wilson loop is the 
same as in general kinematics, although the computation is remarkably simplified with 
respect to general kinematics. Using the simplest of those kinematics, the QMRK of a pair 

3 Note that these expressions are valid for generic values of the parameters a and b. For the limiting cases 
where the parameters approach or 1 some care is needed, e.g., G(0; z) — \nz, whereas lim a _>o In (l — -) 
is divergent. 
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along the ladder, we have performed the first analytic computation of the two-loop six- 
edged Wilson loop in general kinematics. The computation has been done in the Euclidean 
region, where the result is real. Except for the contribution coming from the hard diagram 
with six light-like edges, the result can be expressed in terms of harmonic polylogarithms 
in one conformal cross ratio. In turn, the six-edged hard diagram, which constitutes the 
bulk of the final result, can be written as a linear combination of Goncharov's multiple 
polylogarithms, whose arguments are functions of conformal cross ratios. Finally, the 
remainder function can be expressed as a linear combination of multiple polylogarithms of 
uniform transcendental weight four. Details on the explicit computation of the integrals 
will be presented in a forthcoming publication [47]. 
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